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Fig.  4.  Analysis  of  tissue  sections  from  experimental  shigellosis  in 
rabbits  (A  to  D)  and  spontaneous  human  appendicitis  (E  to  H).  (A) 
immunofluorescence  staining  of  histones  reveals  nuclear  and  extra¬ 
cellular  localization  that  largely  overlaps  with  staining  for  DNA  (C). 
(B)  Staining  with  an  antibody  against  Sh/ge/Za-specific  EPS.  (D)  the 
overlay  indicates  that  numerous  Shigellae  are  closely  associated  to 


fibrous  material  staining  for  histones  and  DNA.  (E)  Staining  for 
neutrophil  elastase  in  an  area  of  neutrophil  exudate  in  human  spon¬ 
taneous  appendicitis  reveals  fibrous  extracellular  material  that  also 
stains  for  histone  (F)  and  DNA  (C).  (H)  Overlay  of  the  images.  The 
images  are  projections  of  confocal  z  stacks  generated  from  sections  of 
5  to  6  p,m  thickness.  Bar,  50  p,m. 
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In  the  Old  World,  some  mosquitoes  in  the  Culex  pipiens  complex  are  excellent 
enzootic  vectors  of  West  Nile  virus,  circulating  the  virus  among  birds,  whereas 
others  bite  mainly  humans  and  other  mammals.  Here  we  show  that,  in  northern 
Europe,  such  forms  differing  in  behavior  and  physiology  have  unique  micro¬ 
satellite  fingerprints  with  no  evidence  of  gene  flow  between  them,  as  would  be 
expected  from  distinct  species.  In  the  United  States,  however,  hybrids  between 
these  forms  are  ubiquitous.  Such  hybrids  between  human-biters  and  bird-biters 
may  be  the  bridge  vectors  contributing  to  the  unprecedented  severity  and  range 
of  the  West  Nile  virus  epidemic  in  North  America. 
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Species  in  the  Culex  pipiens  complex  are 
considered  to  be  the  primary  vectors  of 
West  Nile  virus  (WNV)  in  North  America 
because  they  are  often  the  most  common 
mosquitoes  in  urban  areas  (1),  because  dis¬ 
ease  outbreaks  occur  during  their  peak 
abundance  period  (2),  because  they  are 
competent  laboratory  vectors  of  WNV  (3), 
and  because  field  populations  in  the  United 
States  have  repeatedly  been  found  infected 
with  the  virus  (4,  5).  In  addition,  they  can 
transmit  the  virus  transovarially  (6),  so 


overwintering  mosquitoes  can  serve  as  a 
source  of  WNV  to  initiate  an  infection 
cycle  in  the  spring  (7).  Blood-meal  analysis 
has  revealed  that  Cx.  pipiens  in  the  United 
States  bite  both  humans  (anthropophagy) 
and  birds,  suggesting  they  may  serve  as 
bridge  vectors  of  the  disease  from  birds  to 
humans  (2).  Human  WNV  epidemics  re¬ 
quire  bridge  vectors,  because  humans  and 
other  mammals  do  not  usually  generate 
high  enough  viremia  to  infect  biting  mos¬ 
quitoes  (S).  Although  Cx.  pipiens  has  been 
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implicated  in  recent  deadly  urban  outbreaks 
of  WNV  in  Europe,  these  were  nonrecur¬ 
ring  and  localized  (8).  The  persistence  and 
spread  of  WNV  in  the  United  States  have 
therefore  been  very  surprising,  and  an  in- 
depth  examination  of  different  populations 
of  Cx.  pipiens  worldwide  was  warranted. 

The  Cx.  pipiens  complex  has  been  con¬ 
sidered  “one  of  the  major  outstanding  prob¬ 
lems  in  mosquito  taxonomy”  (9),  because  a 
markedly  divergent  array  of  physiological 
and  behavioral  traits  occurs  without  dis¬ 
tinctive  morphological  differentiation  (1, 
10).  Two  species  are  currently  formally 
recognized  in  the  complex:  Cx.  pipiens  Lin¬ 
naeus  1758  and  Cx.  quinquefasciatus  Say 
1823,  the  northern  and  southern  house  mos- 
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quitoes  that  are  ubiquitous  in  temperate  and 
tropical  regions,  respectively  {11).  Their 
geographical  distributions  overlap  closely 
with  that  of  humans,  who  introduced  them 
into  many  areas  {!). 

Most  specimens  of  the  temperate  spe¬ 
cies,  Cx.  pipiens,  require  a  blood  meal  for 
egg  development  (anautogeny)  and  over¬ 
whelmingly  bite  birds  (72).  In  northern  Eu¬ 
ropean  cities,  however,  some  underground 
populations,  such  as  those  in  the  London 
underground  railway  tunnels  (72),  can  lay  a 
first  batch  of  eggs  without  a  blood  meal 
(autogeny),  although  they  will  bite  humans 
readily  and  ravenously  given  the  opportu¬ 
nity  {13).  These  same  underground  popula¬ 
tions  can  also  breed  in  confined  spaces 
(stenogamy),  whereas  aboveground  Cx. 
pipiens  mate  while  swarming  in  a  large 
breeding  area.  Underground  populations  re¬ 
main  active  throughout  the  year,  whereas 
northern  aboveground  populations  of  Cx. 
pipiens  go  into  torpor  during  the  winter 
(diapause)  (7). 

One  interpretation  of  the  existence  of 
these  two  distinct  strategies  in  Cx.  pipiens 
is  that  rapid  shifts  in  physiological  and 
behavioral  traits  are  the  consequence  of 
local  adaptation  to  underground  environ¬ 
ments  associated  with  human  activity  {14). 
Indeed,  a  study  comparing  underground 
and  aboveground  populations  in  London 
using  allozymes  concluded  that  the  under¬ 
ground  autogenous  populations  were  de¬ 


rived  from  those  aboveground.  However, 
the  authors  found  such  high  differentiation 
and  mating  barriers  between  the  two  phys¬ 
iological  forms  that  they  were  forced  to 
hypothesize  that  a  single  colonization  event 
was  responsible  for  all  the  underground 
populations  in  London  (72). 

An  alternative  interpretation  for  the 
existence  of  two  behavioral  and  physio¬ 
logical  forms  is  that  they  correspond  to  two 
separate  species  (i.e.,  genetically  distinct 
entities).  A  predominantly  northern  spe¬ 
cies,  Cx.  pipiens,  would  be  a  bird- 
dependent  anautogenous  mosquito  that 
diapauses  and  needs  open  spaces  to  mate. 
Its  southern  counterpart,  Cx.  molestus 
[presently  a  synonym  of  Cx.  pipiens  (9, 
73)],  would  not  diapause  and  would  be 
stenogamous.  In  this  scenario,  the  northern 
underground  populations  are  derived  from 
autogenous  Cx.  molestus  of  southern  Eu¬ 
rope  and  northern  Africa,  which  colonized 
northern  countries  as  warm  and  nutrient- 
rich  underground  facilities  were  created 
{14).  In  southern  populations,  both  autog¬ 
enous  and  anautogenous  individuals  co¬ 
exist  (9,  13),  but  autogeny  is  commonly 
associated  with  a  preference  for  mammals 
and  in  particular  anthropophagy  (7,  14). 
Indeed,  anthropophagy  was  the  character 
that  originally  led  to  the  description  of  Cx. 
molestus  Forskal  1775  from  autogenous 
Egyptian  specimens,  although  this  species 
was  later  synonymized  under  Cx.  pipiens 


Fig.  1.  Unrooted  distance  tree, 
based  on  proportion  of  shared 
alleles,  depicting  the  relation¬ 
ships  between  the  33  popula¬ 
tions  used  in  this  study.  Num¬ 
bers  in  parentheses  indicate  the 
number  of  individuals  sampled, 
and  numbers  on  branches  indi¬ 
cate  bootstrap  values.  The  num¬ 
bers  before  each  name  corre¬ 
spond  to  the  numbers  in  Fig.  2. 
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because  of  the  lack  of  diagnostic  morpho¬ 
logical  characters  (13). 

Critically,  these  two  hypotheses  gener¬ 
ate  different  predictions  about  genetic  sim¬ 
ilarity  among  autogenous  and  anautog- 
enous  populations.  If  northern  European 
underground  autogenous  populations 
evolved  from  aboveground  Cx.  pipiens, 
then  they  should  be  genetically  more  close¬ 
ly  related  to  nearby  aboveground  anautog- 
enous  populations  than  to  autogenous 
populations  in  other  geographic  areas.  If, 
however,  underground  autogenous  popula¬ 
tions  belong  to  a  different  species,  then  all 
autogenous  populations  should  constitute  a 
monophyletic  group. 

We  tested  these  hypotheses  using  highly 
polymorphic  microsatellite  markers  (15)  de¬ 
veloped  for  species  in  the  Cx.  pipiens  com¬ 
plex  (16,  17).  We  used  an  array  of  eight 
microsatellite  loci  that  amplify  consistently 
across  all  populations  tested  (CQll,  CQ26, 
CxqGT4,  CxqGTbb,  CxpGT4,  CxpGT9, 
CxpGT12,  and  CxpGT46),  to  fingerprint  au¬ 
togenous  (underground)  and  anautogenous 
(aboveground)  populations  of  Cx.  pipiens 
from  Britain  and  Germany,  as  well  as  other 
aboveground  populations  from  Britain, 
Sweden,  northern  and  southern  France,  It¬ 
aly,  northern  Africa,  the  Middle  East,  Ja¬ 
pan,  Australia,  and  the  United  States  (tables 
SI  to  S2).  Only  autogenous  forms  are 
known  to  have  been  introduced  to  Japan 
and  Australia  (18,  19).  We  also  analyzed 
Cx.  quinquefasciatus  from  southern  locali¬ 
ties  in  the  United  States  as  an  outgroup 
(Fig.  1).  The  similarity  between  autogenous 
mosquitoes  from  different  geographic  areas 
and  their  differentiation  from  sympatric 
or  geographically  close  populations  of 
aboveground  mosquitoes  was  apparent 
from  the  standard  distance  measures  (Fig. 
1),  as  well  as  from  the  allelic  frequency 
distributions  (tables  SI  to  S2).  Northern 
European  aboveground  and  underground 
(autogenous)  populations  formed  two  dis¬ 
tinct  clusters.  Underground  populations 
clustered  with  North  African  and  Middle 
Eastern  populations  as  well  as  with  the 
autogenous  populations  from  Japan  and 


Australia.  Significantly,  the  northern  Euro¬ 
pean  underground  populations  have  a  lower 
average  number  of  alleles  per  locus  than 
the  aboveground  populations  (alleles  ± 
SE  =  2.06  ±  0.20  and  7.64  ±  0.43,  respec¬ 
tively,  P  <  0.01,  averaging  across  both 
loci  and  populations).  Their  most  common 
alleles  are  also  those  that  are  common  in 
African  and  Middle  Eastern  populations 
and  are  often  nonexistent  in  northern 
European  aboveground  populations  (sup¬ 
porting  online  material).  Our  results  there¬ 
fore  reject  the  hypothesis  that  northern 
European  autogenous  mosquitoes  result 
from  local  episodes  of  adaptation  to  under¬ 
ground  environments.  Instead,  they  indi¬ 
cate  that  northern  European  autogenous 
underground  mosquitoes  are  a  separate  ge¬ 
netic  entity,  most  likely  derived  from  a 
southern  species. 

In  the  genetic  distance  analysis,  U.S. 
Cx.  pipiens  populations  clustered  separ¬ 
ately  (Fig.  1).  To  examine  in  more  detail 
the  relationship  between  the  North  Ameri¬ 
can  populations  and  all  remaining  pop¬ 
ulations,  we  performed  a  multilocus  geno¬ 
type  analysis  that  did  not  use  prior 
information  about  the  origins  of  individu¬ 
als.  Instead,  it  combined  all  the  available 
individual  multilocus  genotypes  and  calcu¬ 
lated  the  likelihood  of  a  predetermined 
number  of  clusters  (K),  assuming  that  some 
admixture  (hybridization)  might  take  place 
(i.e.,  we  did  not  assume  genetic  isolation). 
We  obtained  the  most  likely  K  value  from 
the  examination  of  the  range  of  likelihoods 
and  their  distribution  (20)  and  plotted  the 
probability  of  ancestry  from  each  cluster 
for  each  individual  (Fig.  2). 

Using  this  technique,  we  identified  three 
(K  =  3)  most  likely  distinct  genetic  clusters 
in  the  641  specimens  included  in  the  anal¬ 
ysis  (Fig.  2,  the  same  specimens  as  in  Fig. 
1).  These  correspond  to  cluster  A,  which 
includes  the  overwhelming  majority  of 
European  aboveground  populations;  cluster 
B,  which  includes  autogenous  specimens; 
and  cluster  C,  Cx.  quinquefasciatus.  This 
analysis  shows  that  all  U.S.  populations  of 
Cx.  pipiens  include  hybrids,  which  we  de¬ 


fine  as  specimens  with  a  probability  of 
ancestry  equal  to  or  higher  than  0.06  from 
more  than  one  cluster  (more  than  one  color 
in  Fig.  2).  On  average,  more  than  40%  of 
the  specimens  in  each  U.S.  local  population 
had  hybrid  ancestry  from  cluster  A  and 
cluster  B.  No  significant  number  of  hybrids 
was  found  in  northern  European  popula¬ 
tions.  A  small  number  of  hybrids  were 
found  in  two  of  the  six  southern  European 
populations  (10.2%  per  population  on  av¬ 
erage)  (Fig.  2).  Examined  separately,  north¬ 
ern  and  southern  European  populations  did 
not  differ  in  allelic  diversity  (7.64  ±  0.43 
and  7.9  ±  0.46,  respectively),  whereas  U.S. 
Cx.  pipiens  have  a  significantly  lower  av¬ 
erage  number  of  alleles  per  locus  and  per 
population  (5.90  ±  0.30,  P  <  0.01).  The 
U.S.  allelic  depletion  is  consistent  with  the 
idea  that  New  World  populations  of  Cx. 
pipiens  were  introduced  from  the  Old 
World  (1).  We  included  in  the  analyses  a 
population  from  South  Carolina,  a  region 
where  hybridization  between  Cx.  pipiens 
and  Cx.  quinquefasciatus  is  prevalent  (Fig. 
2),  to  show  that  the  uniqueness  of  U.S.  Cx. 
pipiens  populations  is  not  due  to  the  docu¬ 
mented  hybridization  with  Cx.  quinquefas¬ 
ciatus  (21-23).  The  South  Carolina  popu¬ 
lation  clusters  away  from  the  core  U.S.  Cx. 
pipiens  (Fig.  1). 

The  separation  between  the  two  behav¬ 
ioral  and  physiological  forms  of  Cx.  pipi¬ 
ens,  which  for  the  sake  of  brevity  we  will 
herein  refer  to  as  Cx.  pipiens  and  Cx.  mo- 
lestus,  is  probably  quite  recent.  The  micro¬ 
satellite  signatures  of  Cx.  pipiens  and  Cx. 
molestus,  though  distinct,  are  similar  when 
compared  to  Cx.  quinquefasciatus,  for  ex¬ 
ample  (tables  SI  to  S8).  Cx.  molestus  might 
be  a  commensal  form  that  arose  in  as¬ 
sociation  with  the  establishment  of  human 
settlements  after  the  advent  of  agriculture 
(14)  or  by  isolation  during  the  Pleistocene 
glaciations.  Both  scenarios  would  place 
the  separation  at  ~  10,000  years  ago  (24, 
25),  which  could  explain  why  more  slowly 
evolving  allozyme  and  morphological 
traits  do  not  always  distinguish  the  two 
forms  (1,  12,  26,  27). 


Fig.  2.  Results  of  a  Bayesian 
cluster  analysis.  Each  of  the 
641  individuals  included  in  the 
analysis  is  represented  by  a 
thin  vertical  line,  partitioned 
into  three  colored  segments 
that  represent  the  individual's 
probability  of  belonging  to  one 
of  the  three  genetic  clusters 
(yellow,  Cx.  quinquefasciatus-, 
blue,  Cx.  pipiens-,  red,  Cx.  mo¬ 
lestus).  For  this  representation, 
specimens  were  grouped  by  lo¬ 
cation  (bracketed),  and  the  lo- 
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cation  numbers  are  the  same  as  in  Fig.  1.  The  two  long  red  lines  within  southern  France  populations  25  and  27  are  specimens  that  have  a 
probability  of  ancestry  from  Cx.  molestus  of  more  than  0.95  and  are  not  considered  hybrid. 
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What  led  to  the  pervasiveness  of  hybrids 
between  Cx.  pipiens  and  Cx.  moleslus  in 
North  America,  but  not  in  Europe  and  Af¬ 
rica,  still  remains  to  be  determined.  In 
southernmost  Europe,  we  identified  two 
populations  with  a  few  hybrid  individuals, 
as  well  as  populations  with  pure  Cx.  pipiens 
signatures  and  populations  with  a  mix  of 
pure  Cx.  pipiens  and  pure  Cx.  molestus 
signatures  (Fig.  2).  Indeed,  previous 
allozyme-  based  studies  indicated  the  exis¬ 
tence  of  populations  in  Italy  with  a  mix  of 
the  two  forms  (26)  but  a  very  low  rate  of 
hybridization  (1%),  probably  because  of 
their  different  mating  behaviors  (26).  The 
rarity  of  southern  European  hybrids  and  our 
failure  to  find  hybrids  in  northern  Europe 
may  be  due  to  their  low  fitness  and  inability 
to  diapause.  Importantly,  the  introduction  to  the 
United  States  of  separate  populations  of  Cx.  pipi¬ 
ens  and  Cx  molestus  that  later  hybridized,  or  of 
hybrids  from  southern  Europe,  has  led  to  abun¬ 
dant  and  ubiquitous  hybrid  forms  that  survive  the 
rigors  of  northern  winters. 

It  is  now  clear  that  models  derived  from 
the  U.S.  epidemic  of  WNV  (2S)  may  not  be 
applicable  to  Eurasia,  and  vice  versa  (29). 
A  major  factor  in  all  recent  outbreaks  (Ro¬ 
mania  1996,  Russia  1999,  and  United 
States  1999)  is  the  involvement  of  mosqui¬ 
toes  in  the  Cx.  pipiens  complex  as  the 
primary  vectors  {8,  30).  Unlike  European 
Cx.  pipiens,  U.S.  Cx.  pipiens  appears  to  bite 
readily  both  avian  hosts  and  humans  (2, 
31).  Here  we  have  shown  that,  across  the 
northeastern  United  States,  a  large  propor¬ 
tion  of  individuals  are  hybrids  of  human- 
biter  and  bird-biter  forms.  In  combination 
with  susceptible  migrating  birds  and  highly 
concentrated  human  populations  in  U.S. 
cities  and  suburbs,  the  prevalence  of  such 
bridge  vectors  that  readily  transmit  the  vi¬ 
rus  among  and  between  avian  hosts  and 
humans  could  have  created  the  current  ep¬ 
idemic  conditions. 

The  present  study  suggests  that  changes  in 
vectorial  capacity  and  the  creation  of  new 
efficient  vectors  may  occur  with  new  intro¬ 
ductions.  In  particular,  the  arrival  of  hybrid 
American  forms  in  northern  Europe  has  the 
potential  to  radically  change  the  dynamics  of 
WNV  in  Europe. 

References  and  Notes 

1.  E.  B.  Vinagradova,  Culex  pipiens  pipiens  Mosquitoes: 
Taxonomy.  Distribution.  Ecology.  Physiology.  Genet¬ 
ics.  Applied  Importance,  and  Control  (Pensoft,  Mos¬ 
cow,  2000). 

2.  A.  Spielman,  Ann.  N.Y.  Acad.  Sci.  951,  220  (2001). 

3.  M.  J.  Tureli,  M.  L  O'Guinn,  D.  J.  Dohm,  J.  W.  Jones, 
J.  Med.  Entomol.  38,  130  (2001). 

4.  K.  A.  Bernard  et  a/.,  Emerg.  Infect.  Dis.  7,  679 
(2001). 

5.  V.  L.  Kulasekera  et  al.,  Emerg.  Infect.  Dis.  7.  722 
(2001). 

6.  D.  J.  Dohm,  M.  R.  Sardelis,  M.  J.  Tureli,  j.  Med. 
Entomol.  39,  640  (2002). 

7.  R.  S.  Nasci  eta/.,  Emerg.  infect.  Dis.  7,  742  (2001). 


8.  C.  G.  Hayes,  Ann.  N.Y.  Acad.  Sci.  951,  25  (2001). 

9.  R.  E.  Harbach,  C.  Dahl,  G.  B.  White,  Proc.  Entomol. 
Soc.  Wash.  87,  1  (1985). 

10.  T.  Guillemaud,  N.  Pasteur,  F.  Rousset,  Proc.  R.  Soc. 
London  Ser.  B.  264,  245  (1997). 

11.  K.  L  Knight,  Supplement  to  the  Catalog  of  the 
Mosquitoes  of  the  World  (Diptera:  CuUcidae),  T.  S. 
Foundation,  Ed.  (Entomological  Society  of  America, 
College  Park,  MD,  1978),  supplement  to  vol.  6. 

12.  K.  Byrne,  R.  A.  Nichols,  Heredity  82,  7  (1999). 

13.  R.  E.  Harbach,  B.  A.  Harrison,  A.  M.  Gad,  Proc.  Ento¬ 
mol.  Soc.  Wash.  86,  521  (1984). 

14.  P.  F.  Mattingly  et  al..  Trans.  R.  Entomol.  Soc.  Lond. 
102,  331  (1951). 

15.  D.  X.  Zhang,  G.  M.  Hewitt,  Mol.  Ecol.  12,  563 
(2003). 

16.  N.  Keyghobadi,  M.  A.  Matrone,  G.  D.  Ebel,  L.  D. 
Kramer,  D.  M.  Fonseca,  Mol.  Ecol.  Notes  4,  20 
(2004). 

17.  D.  M.  Fonseca,  C.  T.  Atkinson,  R.  C.  Fleischer,  Mol. 
Ecol.  7,  1617  (1998). 

18.  F.  H.  Drummond,  Trans.  R.  Entomol.  Soc.  Lond.  102, 
369  (1951). 

19.  K.  Tanaka,  K.  Mizusawa,  E.  S.  Saugstad,  Contrib.  Am. 
Entomol.  Inst.  16,  1  (1979). 

20.  J.  K.  Pritchard,  M.  Stephens,  P.  Donnelly,  Genetics 
155,  945  (2000). 

21.  A.  R.  Barr,  Am.  J.  Trop.  Med.  Hyg.  6,  153  (1957). 

22.  A.  J.  Cornel  et  al..J.  Med.  Entomol.  40,  36  (2003). 

23.  S.  Urbanelli,  F.  Silvestrini,  W.  K.  Reisen,  E.  De  Vito, 
L.  Bullini,7.  Med.  Entomol.  34,  116  (1997). 

24.  L.  L.  Cavalli-Sforza,  F.  Cavalli-Sforza,  The  Great 
Human  Diasporas:  The  History  of  Diversity  and 
Evolution  (Addison-Wesley,  Reading,  MA,  1995). 

25.  J.  de  Zulueta,  Parassitologia  36,  7  (1994). 

26.  S.  Urbanelli  et  al.,  in  Ecologia.  Atti  I  Congr.  Naz. 


Many  networks  in  nature  share  global  prop¬ 
erties  (I,  2).  Their  degree  sequences  (the 
number  of  edges  per  node)  often  follow  a 
long-tailed  distribution,  in  which  some  nodes 
are  much  more  connected  than  the  average 


Departments  of  Molecular  Cell  Biology,  Physics  of 
Complex  Systems,  and  Computer  Science,  Weizmann 
Institute  of  Science,  Rehovot  76100,  Israel. 

*To  whom  correspondence  should  be  addressed  at 
Department  of  Molecular  Cell  Biology,  Weizmann  In¬ 
stitute  of  Science,  Rehovot  76100,  Israel.  E-mail: 
urialon@weizmann.ac.il 


Soc.  Ital.  Ecol.,  A.  Moroni,  O.  Ravera,  A.  Anelli,  Eds. 
(Zara,  Parma,  Italy,  1981),  pp.  305-316. 

Z7.  C.  Chevillon,  R.  Eritja,  N.  Pasteur,  M.  Raymond,  Genet. 
Res.  66,  147  (1995). 

28.  P.  D.  Crook,  N.  S.  Crowcroft,  D.  W.  Brown,  Commun. 
Dis.  Public  Health  5,  138  (2002). 

29.  Z.  Hubalek,  Viral  Immunol.  13,  415  (2000). 

30.  H.  M.  Savage  et  al.,  Am.  J.  Trop.  Med.  Hyg.  61,  600 
(1999). 

31.  A.  G.  Richards,  Entomol.  Nev/s  52,  211  (1941). 

32.  We  thank  N.  Becker,  C.  Curtis,  M.  Carroll,  G.  Ebel,  C. 
Evans,  M.  Santa  Ana  Gouveia,  L.  Kramer,  G. 
O'Meara,  F.  Noguera,  and  C.  Williams  for  providing 
invaluable  mosquito  samples;  J.  Smith,  M.  Matrone, 
T.  Ganguly,  and  the  DNA  Sequencing  Facility,  Uni¬ 
versity  of  Pennsylvania,  for  technical  assistance; 
and  A.  Bhandoola  and  four  anonymous  reviewers 
for  comments  and  valuable  suggestions  on  an  ear¬ 
lier  version  of  this  manuscript.  Supported  by  a 
National  Research  Council  Associateship  through 
the  Walter  Reed  Army  Institute  of  Research 
(D.M.F.),  by  NIH  grant  nos.  U50/CCU220532  and 
1R01GM063258,  and  by  NSF  grant  no. 
DEB-0083944.  This  material  reflects  the  views  of 
the  authors  and  should  not  be  construed  to  repre¬ 
sent  those  of  the  Department  of  the  Army  or  the 
Department  of  Defense. 

Supporting  Online  Material 

www.sciencemag.org/cgi/content/full/303/5663/1535/ 

DCl 

Materials  and  Methods 

Tables  SI  to  S8 

References  and  Notes 

2  December  2003;  accepted  16  January  2004 


(i).  In  addition,  natural  networks  often  show 
the  small-world  property  of  short  paths  be¬ 
tween  nodes  and  highly  clustered  connections 
(1,  2,  4).  Despite  these  global  similarities, 
networks  from  different  fields  can  have  very 
different  local  structure  (J).  It  was  recently 
found  that  networks  display  certain  patterns, 
termed  “network  motifs,”  at  much  higher  fre¬ 
quency  than  expected  in  randomized  net¬ 
works  (6,  7).  In  biological  networks,  these 
motifs  were  suggested  to  be  recurring  circuit 
elements  that  carry  out  key  information¬ 
processing  tasks  (6,  8-10). 


Superfamilies  of  Evolved  and 
Designed  Networks 

Ron  Milo,  Shalev  Itzkovitz,  Nadav  Kashtan,  Reuven  Levitt, 
Shai  Shen-Orr,  Inbal  Ayzenshtat,  Michal  Sheffer,  Uri  ALon* * 

Complex  biological,  technological,  and  sociological  networks  can  be  of  very 
different  sizes  and  connectivities,  making  it  difficult  to  compare  their  struc¬ 
tures.  Here  we  present  an  approach  to  systematically  study  similarity  in  the 
local  structure  of  networks,  based  on  the  significance  profile  (SP)  of  small 
subgraphs  in  the  network  compared  to  randomized  networks.  We  find 
several  superfamilies  of  previously  unrelated  networks  with  very  similar  SPs. 
One  superfamily,  including  transcription  networks  of  microorganisms,  rep¬ 
resents  "rate-limited"  information-processing  networks  strongly  con¬ 
strained  by  the  response  time  of  their  components.  A  distinct  superfamily 
includes  protein  signaling,  developmental  genetic  networks,  and  neuronal 
wiring.  Additional  superfamilies  include  power  grids,  protein-structure  net¬ 
works  and  geometric  networks.  World  Wide  Web  links  and  social  networks, 
and  word-adjacency  networks  from  different  languages. 
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